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Abstract 

The ability to accurately control a quantum system is a fundamental requirement in many areas 
of modern science [lj such as quantum information processing [2] and the coherent manipulation of 
molecular systems. It is usually necessary to realize these quantum manipulations in the shortest 
possible time in order to minimize decoherence, and with a large stability against fluctuations of 
the control parameters. While optimizing a protocol for speed leads to a natural lower bound in 
the form of the quantum speed limit rooted in the Heisenberg uncertainty principle [3H2] , stability 
against parameter variations typically requires adiabatic following of the system. The ultimate 
goal in quantum control is to prepare a desired state with 100% fidelity. Here we experimentally 
implement optimal control schemes that achieve nearly perfect fidelity for a two-level quantum 
system realized with Bose-Einstein condensates in optical lattices [7j. By suitably tailoring the 
time-dependence of the system's parameters, we transform an initial quantum state into a desired 
final state through a short-cut protocol reaching the maximum speed compatible with the laws 
of quantum mechanics. In the opposite limit we implement the recently proposed transitionless 
superadiabatic protocols [8HTT]. in which the system perfectly follows the instantaneous adiabatic 
ground state. We demonstrate that superadiabatic protocols are extremely robust against param- 
eter variations, making them useful for practical applications 
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The need to control the dynamics of a quantum system is common to many different 
areas of science, ranging from the coherent manipulation of molecular systems [12J to high 
precision measurements [13]. Quantum control may aim at reaching a given target state, as 
in the cooling of atomic ensembles [14] and nano-mechanical oscillators [15j, or at tracking 
the instantaneous ground state of a system during its evolution as in adiabatic quantum 
computation [16]. For each of these tasks an optimum strategy can be designed. In most 
applications the requirements on the fidelity of the final state is stringent, and in principle 
one aims at perfect fidelity. 

Here we investigate high-fidelity quantum control protocols for the 'simplest non-simple 
quantum problem' [17J, i.e. the evolution of a two-level system in a time T as illustrated in 
Figure 1 (energies and times are expressed in the natural units of energy and time of our 
physical system, see Methods). The two diabatic levels |0) and |1) with energies Y{r) = 
±4(r - §), where r = t/T ranges from to 1, are coupled through a coupling parameter uj 
giving rise to the adiabatic levels |^ e (r)). While there are infinitely many paths in Hilbert 
space connecting an initial quantum state 1^^) (for r = 0) with a final target state |\P/m) 
(for r = 1) such that the final fidelity J 7 ^ = \(^ f in \ipg(r = 1))| 2 = 1, in this work we 
concentrate on two special classes. We consider those paths that minimize the time T of the 
transformation and hence reach the quantum speed limit, and those paths ensuring a perfect 
following of the adiabatic ground state \^ 9 {t)). We will call the latter 'superadiabatic' paths 
with reference to Berry's work on transition histories and the superadiabatic basis [8]. 

We realize an effective two-level system using Bose-Einstein condensates (BECs) in an 
accelerated optical lattice (see Methods) . Under appropriate conditions the wavefunction of 
the BEC in the periodic potential of the optical lattices can be approximated by considering 
only the two lowest energy bands [23]. Inside the first Brillouin zone the time-dependent 
Hamiltonian of the system can then be written in terms of the Pauli matrices Oi as 

H = T(r)a z + uj(r)a x . (1) 

The system is initially prepared in the state l^m), and the target is to reach |^fi n ) after an 
evolution of duration T. The initial and final values of the parameters are chosen to be on 
opposite sides of the energy anticrossing (at r = 0.5) (see Fig. 1). 

We begin by considering a protocol that takes the system from |^^) to |^fi n ) with 100% 
fidelity in the shortest possible time T min . By analogy with the equivalent classical case this 
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kind of protocol has been called the 'quantum brachistochrone' P, OS]. If we only impose 
the constraint that uj is constant (otherwise T min as uj — >> oo), we find that the protocol 
shown in Figure 2c minimizes T. This 'composite pulse' protocol (see Methods), in close 
analogy to composite pulses in NMR [19] , represents half a Rabi oscillation with frequency 
uj at r = 0.5, preceded and followed by two short pulses (in theory delta-functions) with a 
pulse area of 7r/4 (see also [20]). The experimentally measured minimum times (Fig. 2d) 
for reaching the target state with fidelity Tj% n ~ 1 approach the quantum speed limit, also 
known as the Fleming or Bhattacharyya bound, given by 

_ arccos|(^ /in |^)| 

^qs — • \£) 

UJ 

One can assess the performance of the composite pulse protocol by comparing it to the 
paradigmatic Landau-Zener (LZ) protocol and to the locally adiabatic protocol proposed by 
Roland and Cerf in the context of adiabatic quantum computation [21] . In the LZ problem 
with constant uj and linearly varying V{r) a system initially prepared in the adiabatic ground 
state undergoes tunneling to the excited state with a finite probability, leading to a fidelity 
J 7 fin = exp(— 2E ^^) < 1. Stricter boundary conditions are imposed in the approach taken by 
Roland and Cerf. On top of the condition that uj be constant they demand local adiabatic 
following to within some small deviation e, i.e. during the entire protocol •F(r) > 1 — e 2 (see 
Methods). In both the LZ and the Roland and Cerf protocols the time to achieve perfect 
adiabaticity, i.e. Tji n = 1, diverges and hence we measure the time to reach Tji n = 0.9 as 
a function of uj instead (shown as a dotted line in Fig. 2e). The results of the comparison 
are shown in Fig. 2d. It is evident that while the LZ protocol is more than an order of 
magnitude slower than the composite pulse protocol for small uj ) the Roland-Cerf protocol 
reaches Tfi n = 0.9 in a time that is only 10 — 50% above the lower limit given by the 
Fleming-Bhattacharyya bound. 

At the opposite extreme of the quantum control spectrum, rather than minimizing the 
total time one can maximize the adiabaticity during the protocol. While the composite pulse 
protocol described above was obtained through optimization given a constraint on uj ) it is 
possible to analytically calculate protocols that ensure T{r) = 1 during the entire evolution. 
The reasoning behind such a transitionless super- adiabatic (or counter- adiabatic [lQl ITT] ) 
protocol is that for a given time-varying Hamiltonian it is always possible to construct 
an auxiliary Hamiltonian % s that cancels the non-adiabatic part of the evolution under 
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T-L alone. It thus ensures transit ionless adiabatic following such that the system evolving 
under T~L + T~L S always remains in the instantaneous adiabatic ground state of T-L with 100% 
probability [El [TO] for a finite duration of the protocol. In general T~L S can be written as 

n s (t) = ihJ2\d t n(t))(n(t)\, (3) 

n 

where \n(t)) are the eigenstates of the original Hamiltonian T~L. For a two-level system of 
the form (1) one finds that 

«.<«)- W 

where (f) = arctan^y. This means that in order to make a given evolution of a two-level 
system perfectly adiabatic one needs to add an interaction term corresponding to a a y Pauli 
matrix [22]. In practice, T~L S can be implemented by introducing an additional interaction 
into the system, e.g. through an extra laser or microwave field. In the case of atoms in 
an optical lattice considered in this work, the additional Hamiltonian can be realized by 
adding a second optical lattice shifted with respect to the first one by is the lattice 

spacing). It can be shown, however, that the effect of this additional field can also be 
achieved through an appropriate transformation r V and uj — >> uj 1 (see Methods), so that 
no extra field is necessary. This result, which is independent of the physical system under 
consideration, means that the resulting protocol is intrinsically more stable as there will be 
no problems associated, e.g, with phase fluctuations between the fields [TUl [TT] . 

For the standard LZ protocol with fixed uj and linearly varying Y{r) the general shape of 
the required transformation is shown in Figure 3a (for details see Methods), together with 
the result of an experiment in which for different sweep durations T the fidelity Tj% n was 
measured both for the linear LZ protocol and for the superadiabatic protocol. We find that 
in the superadiabatic case Tji n ^ 0.98 for all T. Furthermore, a time-resolved measurement 
of T{t) during the sweep (Fig. 3d) shows that the system stays in the ground adiabatic 
state at all times. 

In principle, the superadiabatic transformation of T and uj can be calculated for any sweep 
protocol. There are, however, special cases that are of particular interest. For example, a 
protocol for which the correction in T vanishes, i.e. for which r' = T (except at the beginning 
and the end of the sweep, see Methods), can be expected to be more robust to variations in its 
parameters. This is the case for the superadiabatic 'tangent protocol' (see Methods) shown 
schematically in Figure 3b. As in the case of the linear LZ protocol with superadiabatic 



corrections, in the tangent protocol the systems remains in the adiabatic ground state to 
within 1% throughout the entire protocol. Since the noise in our imaging system does not 
allow us reliably to measure fidelities Tj% n ^ 0.98, we performed an experiment in which 
the tangent sweep was repeated four times, giving an overall probability J 7 ^ « 0.94, from 
which we deduce that, indeed, in a single sweep Tji n ~ 0.99. This figure is compatible with 
a fidelity of 100% as we estimate the non-adiabaticity of the preparation and measurement 
protocols to be on the order of 1%. 

In order to test the sensitivity of the superadiabatic protocols to a (simulated) variation 
in the control parameters, we varied both uj and T around the optimum value and measured 
the fidelity Tji n in each case. The results are summarized in Figure 4a, which shows clearly 
that the tangent sweep is extremely robust with respect to an increase in T or cj, with 
J 7 fin ^0.99 for increases up to 100%. 

Finally, we compare the speed of the superadiabatic tangent protocol with the composite 
pulse protocol as a function of uj' (for the composite pulse protocol uj' = uj being constant). 
Solving Eqn. 10 for T gives a total time for the superadiabatic tangent protocol that depends 
on both uj and uj' (with uj' > uj). It is, therefore, possible to minimize T for a given value of 
uj 1 by choosing an appropriate value for uj. The result of this minimization is shown in Fig. 
4b. Surprisingly, the minimum value of T as a function of uj' for the superadiabatic tangent 
protocol lies below the LZ and Roland and Cerf times and is quite close to the quantum 
speed limit, meaning that the 'penalty' in terms of speed for the requirement of perfect 
adiabatic following is less than one might at first expect. In fact, for small values of uj' one 
can formally let uj —> in the expression for uj' (see Eqn. 21 in Methods), giving T = ^ 
which coincides with the quantum speed limit for orthogonal states with K^/ml^m)! = 0. 

In summary, we have explored high-fidelity quantum control protocols for the evolution of 
an artificial two-level quantum system ranging from the speed-limited to the superadiabatic 
regime. The superadiabatic transformations make it possible to readily implement protocols 
ensuring perfect adiabatic following in a variety of existing applications. In practice, of course 
the choice of protocol will depend on the boundary conditions and physical limitations of 
the system under consideration. If both T and uj can be controlled (to within some limits), 
the superadiabatic protocols provide the possibility of state preparation with 100% fidelity 
with high stability against parameter variations. 
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Methods 



A. Bose-Einstein condensates in optical lattices 

We realize an effective two-level system by loading Bose-Einstein condensates into optical 
lattice potentials [7j of the form ^ cos {^l^xjd^ + 0(t)). In the limit of small lattice depths 
Vo ^5 5 E rec (here the recoil energy E Tec = huo rec = 7r 2 h 2 /2M^, and u) rec = 2tt x 3.15 kHz 
defines the natural units of energy huj rec and time l/uj r ec for our system), the lowest energy 
levels are given by the quadratic dispersion relations of free particles with momenta differing 
by 2hkL coupled through a coupling constant uj = Vq/A. Subtracting the quadratic term in 
the momentum [24] leads to the effective Hamiltonian (1), where the time dependence in 
r is now achieved through a variation of the quasimomentum q in the first Brillouin zone. 
Experimentally, q and hence V = ATiuj rec (q — \) (in this Methods section we shall use the 
explicit physical units wherever appropriate) is controlled through the term (j)(t) which can 
be used to accelerate the lattice, leading to an inertial force in the rest frame of the lattice. 
The time dependence of uj is controlled through the power of the lattice laser beams which 
determine Vq. 

The experimental protocols are carried out using techniques previously developed by us 
and described in detail in [23l [24] . Time-resolved measurements with superadiabatic driving 
are performed by applying appropriate jumps of the lattice position, of the quasimomentum 
and/or of the lattice depth, before measuring in the adiabatic basis (now of the original 
Hamiltonian T~L). 

B. Driving protocols 

1. Composite pulse protocol 

The optimality of the time ^ in the LZ setting was derived in Ref. [3j by starting from 
a guess function similar to the Roland and Cerf protocol (see next section), and by running 
a numerical search aimed at determining the minimal value of T (in that optimization the 
adiabatic requirement was not enforced). The optimal pulse T{r) associated with such a 
minimal time corresponds to highly irregular functions which are strongly peaked at the 
beginning and at the end of the protocol while remaining flat and equal to zero for interme- 
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diate times. Assuming that no constraints are imposed on T{r) and on its first derivative 
one can extrapolate the exact (asymptotic) analytical form of such pulses as 



r(r) 



-r forr = 

r M for r e [0,r ] 

for r e [t ,1-t ] (5) 

-r M for r G [1 - T , 1] 

+r for r = 1 , 



where r = —2 as usual, while Ym and ro are, respectively, asymptotically large and small 
quantities which satisfy the condition 

r M T = tt/4 . (6) 

For Y M 3> 1 the evolution described by the pulse ^ corresponds to first applying a fast 
(instantaneous) clock-wise a z pulse around the z-axis to the system which take the point 
on the Bloch sphere lying in the x — z plane and, in the short time To, rotates it into the 
y — z plane. The system then undergoes a rotation at frequency uu around the x— axis for 
a time 1 — 2r , while finally another instantaneous o z pulse rotates the x — z plane back to 
its initial position. 

The exact transfer of \^i n i) to \^fi n ) is achieved by choosing T such that 

m/ , arccos \(^ fi n \^ini)\ , x 

T(l - 2r ) = lx finl ini/l , (7) 

UJ 

which coincides with the r.h.s. of Eq. ^ by letting r —> while keeping ([6]). We note 
here that while the fact that the composite pulse protocol realizes that quantum speed limit 
suggests that this time is optimal, a formal proof for this is still missing. 



2. Roland- Cerf protocol 

In the following we summarize the analysis of Roland and Cerf [21 J adapting it to the case 
of the Hamiltonian % of Eq. (1) under the assumption that the coupling u is kept constant. 
In this protocol the requirement is to keep the evolution adiabatic in each infinitesimal time 
interval. At any instant of the evolution the fidelity of the state with the instantaneous 
ground state is required to be .F(t) = \(^P(T)\ip g (r))\ 2 = 1 — e 2 , from which it follows that 



8 



the total time of the protocol is 



T(e) 



1 1 



(8) 



for a time-dependence of Y{t) of the form 



r(r) 



4aj 2 T(e)(r- 1/2) 



(9) 



^l-16e 2 cj 2 T(e) 2 (r-l/2) 2 ' 



(as shown in Fig. 2b in the main text). 
3. Superadiabatic protocol 

In order to implement superadiabatic (transitionless) driving, we recast the Hamiltonian 
U{t)+U s {t) = Y{r)a z + uo{r)(j x + \^a y in the form W = F \r)a z + uj' \r)a x , eliminating the 
need for an extra potential that realizes the a y -term. The necessary transformations r — >• r 7 
and uj —> uj f can be derived by observing that for the physical system used in this paper, 
i.e. matter waves in a periodic potential, the Pauli matrices a x and a y are simply related by 
the spatial displacement operator U&(5x) with Sx = cIl/A, where is the lattice constant. 
In the basis of the plane waves exp (±ikx) defining the two-state subspace {|0), |1)}, this 
operator can be written as 



This means that the a y -teim in the transitionless driving protocol can be realized by adding 
a second periodic potential shifted by rf L /4 and accelerated in the same way as the first one. 
Since the sum of two periodic potentials of the same periodicity is again a periodic potential 
with a modified phase and amplitude, we can write the combined potential 




(10) 



and hence 



U\(d L /4)a x U d (d L /4). 



(11) 




(12) 



as 




(13) 
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where a = |^ and f3 = arctan From this equation, the transformation rule for uj = Vo/4 
can be read off immediately, while the time derivative of the displacement term 4^ gives 
the correction to the quasimomentum 

4 



q'(t) = q(t) - 



8uo ri 



(14) 



from which the transformation rule for T is calculated using the relation T = 4hu rec (q 
The complete transformation then reads (in rescaled units) 



r' = r- 



1 d 

2dt 



arctan 



/ Tu — uT 
\2u(T 2 + u 2 ), 



7 {Tu - ut) 2 
\ + 4oj 2 (T 2 + oj 2 ) 2 



(15) 



(16) 



where the discontinuities in the expression in square brackets in V at the beginning and at 
the end of the protocol give rise to delta-functions, which can be realized in practice using 
large but finite corrections AI^m for a short duration At such that 



AtAT M = T - 2 arctan (^^y 



(17) 



where the — and + signs refer to the corrections at the beginning and at the end of the 
protocol, respectively. 

The two superadiabatic protocols considered in this paper are the superadiabatic linear 
protocol with 



T'(r) = r(r) - 
u'{t) =ufl 



4(r-|) 



T2[(r-|)2+Ia»2] +l' 



T 2 (8(r-|) 2 + | W 2 ) 2 



(18) 
(19) 



and the superadiabatic tangent protocol for which 

r'(r) = r(r) = cjtan ^2 (r - ^ arctan Q^j^j (20) 
(apart from the delta-functions at the beginning and at the end of the protocol) and 



uj = UJ 



i 



i+ 



arctan(-) 



2^2 



(21) 



(Tu) 2 ■ 

The latter protocol is found by demanding that T'(t) = T{r) and solving the resulting 
differential equation. 
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Figure captions 

FIG. 1: Schematic of a two-level quantum system with crossing energy levels. The bare (diabatic) 
states |0) and |1) are coupled to give the adiabatic states |\E^) and |\I> e ). The adiabatic states have 
an energy gap of 2uj where the diabatic levels cross (at r — 0.5, E = 0). The arrows indicate the 
two extreme optimized protocols discussed in this paper: the 'short cut' composite pulse protocol 
reaching the quantum speed limit (blue line) and the superadiabatic protocol (red line) for which 
the system perfectly follows the instantaneous ground state. 

FIG. 2: Comparison between the Landau Zener (a), Roland and Cerf (b) and composite pulse 
protocols (c). (d) The time needed to achieve a fidelity of 100% for the composite pulse protocol 
(red triangles) and the minimum time to achieve Tf{ n = 0.9 in the Roland-Cerf protocol (grey 
circles) and the LZ protocol (blue squares), (e) Fidelity of the final state as a function of the 
duration for the composite pulse protocol (red triangles), the Roland-Cerf protocol (grey circles) 
and the LZ protocol (blue squares). Note that Tji n — 1 is not reached in the composite pulse 
protocol because of an interaction-induced loss of coherence in the BEC, which also leads to fidelities 
that are slightly below the theoretical predictions for the other protocols. The dashed lines are 
theoretical predictions. All experimental data are for uj — 0.5. 

FIG. 3: Superadiabatic dynamics in a two-level system. (a,b) Original (dashed lines) and supera- 
diabatic (solid lines) protocols for the linear (a) and tangent (b) cases, (c) Fidelity as a function 
of the duration of the protocol for the superadiabatic linear (grey circles) and LZ protocols (blue 
squares). For comparison, the results of a linear LZ sweep in which only uj is transformed (Eq. 8) 
is also shown (empty squares). The dashed lines are theoretical predictions, (d) Fidelity during 
the protocol for the superadiabatic linear (grey circles), superadiabatic tangent (red squares) and 
LZ protocols (blue squares). Clearly, for the superadiabatic protocols .F(t) ~ 1 during the entire 
protocol (to within the experimental error). All experimental data are for uj — 0.55. 
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FIG. 4: Robustness and speed of the superadiabatic tangent protocol, (a) Fidelity of the supera- 
diabatic tangent protocol as a function of the relative deviation of T (red squares) and uj (empty 
squares) from their ideal values uj — 0.5, T — 5.9. The solid and dashed lines are numerical simula- 
tions. Inset: Detail of the graph, highlighting the stability of the protocol for AT/T, Auj/uj > 0. (b) 
Comparison of the duration of superadiabatic tangent protocols with the composite pulse (dashed 
red line) , Roland and Cerf (dashed black line) and Landau Zener protocols (dashed blue line) . The 
solid black line shows the minimum time for the superadiabatic tangent protocol as a function of 
J. 
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